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Abstract. Several models based on the diffusion-limited aggregation (DLA) model were proposed and their 
scaling properties explored by computational and theoretical approaches. In this paper, we consider a new 
extension of the on-lattice DLA model in which the unitary random steps are replaced by random flights 
of fixed length. This procedure reduces the screening for particle penetration present in the original DLA 
model and, consequently, generates new pattern classes. The patterns have DLA-like scaling properties at 
small length of the random fiights. However, as the flight size increases, the patterns are initially round 
and compact but become fractal for sufficiently large clusters. Their radius of gyration and number of 
particles at the cluster surface scale asymptotically as in the original DLA model. The transition between 
compact and fractal patterns is characterized by wavelength selection, and 1/k noise was observed far from 
the transition. 
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\ 1 Introduction 



' Pattern formation is a field of great interest in the Non- 
\ equilibrium Statistical Physics. In special, the diffusion- 
. limited aggregation (DLA) model T is a noteworthy ex- 
' ample in which a very simple algorithm generates com- 
\ plex disorderly patterns. This model was related to several 
' physical and biological applications, such as electrodeposi- 
\ tion viscous fingering 3 , bacterial colonies neurite 
. formation j5j , and, more recently, tumor growth ■ Af- 
' ter two decades, the DLA scaling properties yet remain 
, not completely understood and this model constitutes an 
■ actual and challenger theoretical problem. In the original 
\ model ^ , the simulation begins with an initial seed at the 
. center of the lattice. A particle, represented by a site of the 
' lattice, is released at a random distant point (the launch- 
, ing radius) and performs an on-lattice random walk. If the 
particle visits a neighbor site of the initial seed, it joins ir- 
reversibly to this site. If the distance of the particle to the 
cluster is too large, the particle is excluded and a new one 
is released from the launching circle. Successive new par- 
ticles are released from a new random distant point and 
walk at random until they find the cluster neighborhood. 
This algorithm is inefficient and clusters containing just 
a few thousand particles can be generated with a consid- 
erable computational effort. However, much more efficient 
algorithms can be used in order to grow clusters with more 
than 10® particles for lattices [SIE] and off-lattices pnii| 
DLA models. 



Due its importance as a fundamental model, several 
variants of the DLA model were proposed jjl]. In spe- 
cial, models focusing its screening properties had remark- 
able interest. For example, Meakin investigated a model in 
which the random walks (trajectories with fractal dimen- 
sion 2.0) of the DLA model are replaced by fractal trajec- 
tories (Levy flights or Levy walks) As the fractal di- 
mensions of the trajectories decrease, the particles become 
more penetrating and, consequently, the DLA screening 
effects less intensive. As a consequence, the fractal dimen- 
sions of the clusters increase and their branches become 
denser. Also, drift-diffusion-limited aggregation models, in 
which a particle follows radially biased random walks to- 
ward 1131 and away |14| from the initial seed particle, were 
considered. The former case generates DLA-like clusters 
on shorter length scales and circular-dense patterns, with 
fractal dimension 2.0 on longer length scales as the drift 
probability increases. The later case also generates DLA- 
like clusters on shorter length scales, but eccentric pat- 
terns are generated on longer length scales, with fractal 
dimensions close to 1.0. Moreover, in both cases the pat- 
terns are sparse containing several internal holes. In order 
to generate viscous fingering patterns, several models in- 
troducing surface tension and surface relaxation effects 
in the DLA model were proposed |15II16II17II18| . In these 
models, a sticking probability dependent of the local cur- 
vature for the particles arriving at the surface of the ag- 
gregate is assumed. If the particle does not stick to the 
cluster, it continues diffusing. Again, the DLA screening 
is weakened since the particles do not stick with proba- 
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bility 1. Depending on the surface tension, these models 
generate Eden-hke round and compact patterns '19' on 
shorter length scales that become branched for asymptot- 
ically large clusters. 

In this paper, we present a new extension of the lat- 
tice DLA model in which the unitary random steps are 
replaced by random flights of fixed length. This procedure 
allows the particles to penetrate the screening present in 
DLA clusters and, consequently, generates new pattern 
classes. In Sec. |21 the model and the computational algo- 
rithms are described. In Sec.|31the model results are pre- 
sented and discussed. Finally, some conclusions are drawn 
in Sec. El 

2 Model and methods 

In the modified DLA model presented in this paper, the 
particles perform off-lattice random flights whose step length 
is given by zi = Sa, where a is the lattice constant (as- 
sumed unitary) and 5 > 1 is the unique model parameter. 
Notice that flights substitute walks and, consequently, at 
each step the particle occupies a random position in a 
hypersphere of radius S centered on its initial position. 
Thus, depending on the S value, a particle near to the 
border of the aggregate can reach its inner regions forbid- 
den to particles executing normal random walks. For sake 
of simplicity, the particles of the aggregate lie on a regular 
square lattice. Two conditions determine when the flight- 
ing particles flnish their trajectories becoming part of the 
aggregate. If a particle visits a nearest neighbor site of the 
aggregate, it sticks irreversibly to this site (like in the DLA 
model). But, if the particle reach a site belonging to the 
cluster, it realizes on-lattice random walks of unitary step 
until find an empty site, irreversible sticking on it. Here, 
a surface relaxation, also present in the sticking proba- 
bility dependent models ,15.16, 17..18) . was introduced. In 
order to determine when a particle visits a site, its lattice 
coordinates are defined as the nearest integers associated 
to its real coordinates. The particular case i5 = 1 must 
reproduce the semi-lattice DLA model 20^. 

To optimize the algorithm efficiency some procedures 
were adopted. The random trajectories started at a launch- 
ing circle of radius Rstart = i?o + 55 centered on the initial 
seed. Here, i?o is the maximum distance of a site of the 
aggregate from the center of the lattice and the quan- 
tity 5(5 warranties that the particle realizes some steps 
before it reaches the aggregate. This procedure is justified 
since a random walker released far outside the launch- 
ing circle will first intercept it at a random position with 
equal probability ^J. A kiUing radius Rkui = 100(5 i?o, 
proportional to the flight length, was defined. If a random 
walker cross the circle of radius Rkui centered on the initial 
seed it is discarded. Also, we used a standard method to 
speed up the simulations |1J,J , which consists in allow the 
particles outside the launching circle take long steps ^out 
(longer than 5). However, these long steps cannot bring up 
the particle to inside the launching circle. Thus, the step 
length was chosen as ^ont = "ca&'y-iR — Rstart ~ 5(5, 6), where 
R is the distance of the random walker from the center of 



the lattice. Following the last definition, one can see that 
Syout > S only for R > Rstart + 5(5. In addition, even inside 
the cluster delimiting circle exit large empty regions. So, 
we adapted a algorithm used by Ball and Brady jH|, in 
which the random walker jump a distance if the re- 
gion delimited by the circle of radius centered on the 
walker does not contain occupied sites nor sites neighbor- 
ing the cluster. Thus, the radius must define a region 
containing only sites from which the random walker can- 
not reach the aggregate by a single step. Therefore, 
was determined by the following procedure. The border of 
a sub-lattice of linear size I = max(40, i?o/10) centered on 
the random walker is visited. If none part of the aggregate 
was found, then the sub-lattice is empty and = //2 is 
adopted. Otherwise, if any part of aggregate was found, / 
is reduced by a half and the whole procedure is repeated 
until found an empty sub-lattice or its size reaches the 
value / < 56. In former case a ^ = is assumed and in 
the latter a = (5 is adopted. Fig. ^ shows a scheme of 
the methods described in this paragraph. 




Fig. 1. Optimized DLA model. All the radius defined in the al- 
gorithm are shown: the delimiting aggregate (-Ro), the launch- 
ing (Rstart) and sinking (Rkui) circles, as well as the circles 
that illustrate the long flights used when the particles are dis- 
tant from the aggregate whether inside (Cin) or outside {^out) 
the launching circle. The center of the lattice is depicted in 
black and the other sites of the aggregate are represented by 
white squares. Also, two trajectories, one toward the aggregate 
and the other in which the particle reaches the killing radius, 
are drawn. 



As the step length S increases the patterns become 
denser and, consequently, the random walks over the ag- 
gregate are very computational time consuming. Thus, for 
(5 > 20, the random walkers over the occupied regions were 
also allowed to realize long steps if they are far away from 
the border. 
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3 Results and Discussions 

The DLA model with screening breakdown was simulated 
in square lattices containing up to 12000 x 12000 sites 
and flight lengths ranging in the interval 1 < S < 100. 
The maximum number of particles in the aggregates var- 
ied from N ~ 10^, for smaller S values, to ~ 10^, for 
the larger ones. In Fig. |21 growth patterns generated by 
the model with distinct 6 values are shown. As the flight 
length increases from S — 1, the branches of the clusters 
become denser and for larger S values {S = 40) compact 
patterns, with a nearly regular shape, emerge. In order to 
quantify the geometry of the clusters, their radius of gyra- 
tion Rg and the number of particles on the cluster border 
S were evaluated. These quantities scale with the number 
of particles as 

Rg ^ A^" (1) 

and 

S - A^'^. (2) 

Also, the fractal dimension is given hy df — 1/v. As ex- 
pected, the particular case 5—1 reproduces the semi- 
lattice DLA model. The fractal dimension obtained from 
the slope c? log A^/d log i?g is df = 1.69 ± 0.03, a value 
smaller than df = 1.715 ± 0.004 found for large off-lattice 
DLA simulations [Hj. However, it is a well-understood fact 
that anisotropy effects of the square lattice reduce the ef- 
fective fractal dimension of larger clusters Indeed, 
simulations using a recent method proposed by Bogoy- 
avlenskiy [2] to avoid the lattice anisotropy for DLA 
clusters provide a fractal dimension df = 1.716 ± 0.022 
for 5—1 corroborating our previous claim. 

New patterns emerge due to the weakness of the DLA 
screening resultant from the increasing of the flight lengths. 
This phenomenon can be qualitatively understood through 
an analogy with the sticking probability DLA models. In 
such models there is a correlation length A that increases 
as the sticking probabilities decreases. Thus, the pat- 
terns are compact on smaller length scales (£ ^ A) but 
have DLA-like scaling properties at longer length scales 
{i ^ A) ^J. In order to use this argument for the DLA 
model with screening breakdown, we consider that the ef- 
fective penetration depth of the particles in the sticking 
probability DLA models is proportional to ^JljR,. In- 
deed, IjPa is the mean number of steps realized by the 
particles in the nearness of the cluster border before they 
attach to the cluster and the mean square displacement of 
the random walker is proportional to the number of steps. 
This depth can be associated with the size of the flights 
5 and, consequently, the previous arguments deduced for 
the sticking probability problem can readily be extended 
to the present model. 

The above mentioned length scales are neatly observed 
in Fig. 121 in which the surface sites of the growing clusters 
for distinct numbers of particles are drawn. As one can 
see, the more inner boundaries of Figs. |2Ia) and (b) are 
initially Eden-like, exhibiting a round and rough border. 
As new particles are added to the clusters, their inter- 
faces become progressively regular and, at a given stage, 
quasi-periodic finger-like structures are build up. At this 




Fig. 2. Growth patterns generated by the DLA model with 
screening breakdown. These clusters were generated in square 
lattices containing 2000 x 2000 sites. The aggregate with fewer 
particles (5 — 1) contains 1.5 x 10^ particles whereas the one 
with the larger number (5 = 60) contains 2.8 x lO'^ particles. 



stage, the patterns are featured by a characteristic wave- 
length spontaneously selected. However, this is a transient 
behavior because the "sprouts" consecutively split gener- 
ating fractal structures at long times as illustrated in Fig. 
IHl^c). In this figure, the sequence of the simulation shown 
in Fig.|3Ib) after three doubhngs is drawn. Thus, any com- 
pact pattern becomes fractal for sufficiently large clusters. 
This behavior is analogue to that found in the DLA mod- 
els with sticking probabilities |15lll6j , but the layer mech- 
anisms are different from those involved in such models. 
There, the local curvature of the surface, that is explicitly 
included in those models, seems to develop a essential rule, 
while in the present model only the screening breakdown 
is taken in to account. 

The radius of gyration Rg and the number of periph- 
eral particles S are able to detect the crossover between 
Eden-like and DLA-like patterns suggested by Fig. |3| In 
Fig. Efa) , plots of 5 as a function of the number of par- 
ticles A^ for distinct fiight length sizes 5 are shown. For 
small 5 values an unique scaling law, in which 5* ~ A^, is 
observed indicating that these patterns scale as the DLA 
model. As 5 is increased, the curves exhibit two distinct 
slopes separated by a neat crossover. Below a characteris- 
tic number of particles A^x , the clusters scale as the Eden 
model, i. e., S" N^l"^ ^ and scale as the DLA model above 
this value. The simulations show that for 5 > 8.0 A^x in- 
creases with 5 as A^x ~ 5°^, with a = 3. In the inset of 
Fig. 2fa) , the curve A^x against 5 and the correspondent 
power law fitting are shown. It is evident the quality of 
the fitting. For 5 < 8.0 the determination of the crossover 
becomes difficult. An identical crossover was observed for 
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Thus, we have that zi = aa — 7 = 3/2 and Z2 = av — 7 
0.249. 



Fig. 3. Cluster boundaries along the growth of the aggregate. 
The central boundary represents a cluster containing 10* parti- 
cles and the subsequent ones correspond to patterns after each 
cluster size doubling. Patterns with 5 x 10® particles for (a) 
5 = 20 and (b) 5 = 40 are shown. In (c) the continuation of 
the simulation for 5 = 40 up to 3 x 10^ particles is shown. The 
figure (c) was reduced 50% if compared with the figures (a) 
and (b). 



the radius of gyration. Moreover, the number of peripheral 
sites and the radius of gyration do not depend on the 5 val- 
ues for N < Nx , as indicated in Fig.Efa). However, above 
the crossover, both quantities increase as 6 decreases. In- 
deed, we have that S ^ i5~^^ and Rg S~^^. In figures 
\Mh) and He) the curves S/S'^ and Rg/S^ against iV/J", 
respectively, are plotted. The 7 exponent was varied in 
search for a very satisfactory collapse, and the best value 
is 7 « 3/2. Using that the collapsed curves scale as in the 
Eden model for N ^ Nx and as in the DLA model for 
N ^ Nx , the following scaling relationships are defined 



S^N^f 




(3) 
(4) 



where a and v are the exponents defined in Eqs. ^ and 
10) for the original DLA model, i. e., cr = 1 and u « 0.583. 
The scaling function f{x) is defined as 




ni| — I I I n i i ij — I I I — I I 1 1 i nn — I 1 1 iiiii| — I I 1 1 

lE+1 lE+J lE+3 lE-M lE+5 lE+4 lE+7 

N 




const, if X <C 1 
X if X ^ 1 



(5) 



Nig 



Fig. 4. Radius of gyration Rg and the number of peripheral 
sites 5* as a function of the total number of particles TV in the 
cluster. In (a) the curves S xN are shown for different 5 values. 
In (b) and (c) the collapsed curves of 5* and Rg, respectively, 
are shown for S ranging from S = 9.0 up to 5 = 100. 



One remarkable feature of the present model is the 
mode selection exhibited by the patterns at the compact 
to fractal transition illustrated in Figs. and 0] In order 
to analyze this feature, firstly we map the cluster contour 
in a 1 -|- 1 profile, in which the heights represent the dis- 
tance of a point in the contour from the center-of-mass of 
the pattern. Secondly, we calculated the Fourier spectra 
S{k) of the profile by using the standard FFT method. 
Here, k is the wavenumber and l/k the correspondent 
wavelength. In Fig.^l the mapped profiles of the contours 
drawn in Fig. |31 and the correspondent Fourier spectra 
are shown. For the profile correspondent to the pattern 
before the crossover (Figs. |SJa) andlH^d)) no mode is se- 
lected and the observed low amplitude wavenumbers in- 
dicate the onset of the instabilities that will lead to the 
fingering shown in Fig.lSJb). The Fourier spectrum of the 
cluster border when TV « Nx is characterized by the selec- 
tion of a principal mode and additional lower-amplitude 
modes as indicated in Fig.[SJe). As the cluster increases, 
the selected wavenumber shifts to the left, indicating a 
larger wavelength selection, and additional modes emerge 
in the Fourier spectrum. This results arc in qualitative 
agreement with fluid-fluid displacement experiments for 
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Fig. 5. Profiles of (a) tlie last but two and (b) the last contours shown in Fig. OJb). (c) Profile of the last contour of Fig. 
Etc). The correspondent Fourier spectra are shown at right ((d)-(f)). The profiles are drawn in the same vertical scale whereas 
the horizontal ones depend on the contour. In turn, the Fourier spectra were plotted by using the same horizontal and vertical 
scales. 



radial fingering patterns [22) and recent experiments of 
grain-grain displacement in a Hele-Shaw cell |2Hj . 

In the DLA-like scaling regime {N >> A'^x) a wide 
range of wavelengths are selected as suggested in Fig.|SIf). 
This indicates the presence of self-affine profiles character- 
ized by l/fc^ noise, i. e., all wavelengths are present in the 
Fourier spectra. Thus, the power spectrum of the profile 
scales as 

\S{k)\'^k~^, (6) 

where ^ = 2H + 1 establishes the relationship between ( 
and the Hurst exponent H JT]. Since the power spectra 
are noisy and the wavenumbers are, due the FFT regres- 
sion, evenly spaced, it is hard to determine the power law 
exponent. Therefore, we were led to consider the band- 
integrated regression |24| . in which the power spectrum 
estimated by FFT is integrated over the wavenumber in- 
tervals [fco,/3fco], [(3ko, P'^ko], [/3^fco, /3^fco], and so on. So, 
we computed the quantity 



I{P) = 



\S{k)\ 



(7) 



where kp = /^p/cq and t? ee 1 - C = H/2. In Fig. the 
integrated power spectrum for simulations with (5 = 10 is 
shown. The 77 exponent varies during the growth of the ag- 
gregate from the value t] = 2 {H — 1), that characterizes 
a flat profile, to the value 77 w 1.35 (i7 « 0.68), indicating 
a positively persistent profile, far from the crossover. Our 
simulations were not able to determine the asymptotic 



7] values for all distinct 5 values studied due to compu- 
tational limitations (the maximum lattice size simulated 
was 12000 X 12000). However, simulations for 5 G [8,18] 
provide approximately the same asymptotic rj value and, 
consequently, we conjecture that this exponent is inde- 
pendent of 6. It is important to notice that for 5 < 8 
the pattern contours cannot be mapped the in self-affine 
profiles because very thin branches are present and the 
sequential mapping algorithm used becomes undefined. 




Fig. 6. Typical integrated power spectrum for profiles gen- 
erated from the pattern contours. The flight length used was 
5 = 10 range of and the fit is excludes very long and very short 
wavelengths. 
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An important feature of the DLA model with screening 
breakdown is the formation of isotropic patterns for large 
6 values. This is in marked contrast with the anisotropic 
patterns generated for smaller S values, a central feature 
of on-lattice DLA models m| (see Figs. Inland |3Jl. In or- 
der to quantify the cluster anisotropy, we calculated the 
fourth circular harmonic (cos 40), where is the angu- 
lar position of the particle measured from the initial seed 
position. The averages were done over all peripheral par- 
ticles of up to 100 samples. This procedure excludes the 
effects of internal compact regions. For the original DLA 
model in a square lattices, the fourth harmonic continu- 
ously increases as new particles are added to the cluster 

whereas a null value is found for isotropic patterns. 
In Fig. (cos 4:9) is plotted as a function of the number 
of particles. For smaller 6 values {S < 16), the clusters 
are initially isotropic, i. e., (cos 40} « 0, but an effective 
increasing anisotropy emerges during the cluster growth. 
Such finds agree with the results obtained for the original 
DLA model jQ. Nevertheless, concerning our simulations, 
for larger 6 values {S > 18) none effective anisotropy were 
measured. One can argue that this result may be a con- 
sequence of our computational limitations, i. .e, it did not 
allow us to generate sufficiently large clusters for which a 
measurable anisotropy is present. However, we conjecture 
that these patterns are isotropic since the characteristic 
size for which the lattice anisotropy emerges grows con- 
tinuously up to 5 « 16 and suddenly diverges for S > 18. 
Moreover, the clusters scale as the off-lattice DLA model 
in the asymptotic limit, as shown in Fig. 01 reinforcing the 
idea that lattice effects are not present. The physical ori- 
gin of anisotropy break is related to the relaxation process 
of the particles when they reach the inner regions of the 
cluster. 



4 Conclusions 

In the present work, it was studied an extension of the 
diffusion-limited-aggregation (DLA) model, in which the 
screening effects were directly considered. The unitary steps 
used in the random walks of the original DLA model were 
replaced by random flights of fixed length S. An increasing 
6 leads to a decreasing screening for the penetration of the 
particles in the inner regions of the clusters. 

The model generates DLA-like structures for small S 
values. However, for larger S values the patterns have the 
Eden model scaling properties on smaller length scales 
but scale as the DLA model on longer length scales. The 
compact-to- fractal transition was characterized by the ra- 
dius of gyration Rg and the number of particles on the 
border of the cluster S. Both quantities exhibit two neat 
scaling regimes separated by a characteristic number of 
particles , i. e., the patterns scale as the Eden model for 
iV < iVx but scale as the DLA model for > TVx • From 
a scaling analysis we found that ^ (5^", S ^ (5"^^ and 
Rg ~ (5~^2 with a « 3.00, zi w 1.50 and « 0.249. 

When N « , the patterns exhibit a quasi-periodic 
finger-like structure characterized by the selection of a 
principal and other additional lower-amplitude modes in 
the correspondent Fourier spectra of the cluster bound- 
ary. In turn, in the asymptotic DLA-like scaling regime 
all wavelengths are present and, consequently, a self-affine 
boundary, with Hurst exponent H « 0.68 apparently in- 
dependent on (5, was observed. Moreover, the simulations 
strongly suggest that the lattice effects are not present for 
large S values. 
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Fig. 7. Mean values of the fourth circular harmonic for distinct 
6 values. The number of samples used on the averages varies 
from 100 for 5 = 4 to 30 for 5 = 24. 
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